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^ , We present comprehensive results for the spectrum and decay con- 

stants of hadrons containing a single b quark. The heavy quark is simu- 
lated using an 0(1/M) NRQCD action and the hght quark using the 0{a) 
tadpole-improved clover action on gauge configurations containing two de- 
generate flavours of sea quarks at P^f^'^ = 5.6 provided by the HEM- 
CGC collaboration. We present detailed results for the lo-wer lying S and 
P wave B meson states and the baryon. We find broad agreement 
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with experiment. In addition, we present results for the pseudoscalar and, 
for the first time, the vector decay constants fully consistent to 0{a/M): 
fB = 186(5)(stat)(19)(pert)(9)(disc)(13)(NRQCD)(+50)(a-i) MeV, f% = 
181(6) (stat)(18)(pert)(9)(disc)(13)(NRQCD)(+55)(a-i) MeV and /bJ/b = 
1.14(2) {stat)(—2)(K,s). We present an investigation of sea quark effects in the 
B spectrum and decay constants. We compare our results with those from 
similar quenched simulations at f3^f^^ = 6.0. For the spectrum, the quenched 
results reproduce the experimental spectrum well and there is no significant 
difference between the quenched and Uf = 2 results. For the decay constants, 
our results suggest that sea quark effects may be large. 

I. INTRODUCTION 

Hadrons containing a single b quark hold the key to many important questions facing 
particle physics. In particular, the weak decays of these particles are being studied to look 
for inconsistencies in the Standard Model and indications of new physics beyond it. Lattice 
calculations have a central and fundamental role to play in this pursuit. Not only is lattice 
theory a first principles approach but also offers the most reliable way to calculate the 
masses of heavy-light hadrons and the low-energy QCD factors which are needed to extract 
the electroweak physics from experiment. An introduction to B physics and the theoretical 
advances in this field, for example Heavy Quark Effective Theory (HQET), as well as the 
contribution made from lattice theory, can be found in reference 

The aim of a lattice calculation of B meson phenomenology is to produce reliable predic- 
tions where the systematic errors, such as finite lattice spacing, finite volume and quenching, 
are understood and under control. Our approach is to systematically improve the actions 
used for the heavy and light quark and to use lattices with physical volumes large enough 
to accommodate at least the lower lying B mesons. In particular, we simulate the b quark 
on the lattice using NRQCD: a cut-off of the order of the heavy mass is imposed and we 
cannot extrapolate to the continuum. This requires the systematic errors to be reduced to 
the order of the statistical errors at finite p. In this paper, we argue that the uncertainty 
arising from finite a, volume, and the truncation of the NRQCD series are under control 
and it is now reasonable to investigate the effects of quenching. As a first step towards 
predicting i?-meson properties in full QCD, we perform a simulation of the B spectrum and 
decay constants including two degenerate flavours of dynamical quarks. We then make a 
comparison with our previous NRQCD calculations on quenched configurations and 
study the sensitivity of various quantities to the presence of sea quarks. 

The paper is organised as follows: in the next section we describe the details of the 
simulation. Each choice of quark action and simulation parameter leads to an associated 
systematic error and the corresponding effect on predictions for the B spectrum and decay 
constants are discussed in section [IT^ Results for the spectrum are presented in section 
and compared with experiment. We then study the effect of sea quarks on the B spectrum 
and Af, in section [IV A . 

The next section deals with the pseudoscalar and vector decay constants. The lattice 
operators required to compute the decay constants to 0(1/M) and the corresponding 1-loop 
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perturbative matching factors are introduced and discussed. In addition, the expectations 
from Heavy Quark Symmetry (HQS) for the heavy quark mass dependence of these quantities 
is outhned. Our results for fs, fssi fs* and fs* are presented in section |V A| and their 
dependence on the heavy quark mass is investigated. Finally, we present a detailed study 
of sea quarks effects in the decay constants. 



II. SIMULATION DETAILS 

The simulations were performed using 100 16^ x 32 gauge configurations at f5 = 5.6 with 
two flavours of staggered dynamical sea quarks with a bare quark mass of am^ea = 0.01, 
which roughly corresponds to mT^/mp = 0.525. These configurations were generously made 
available by the HEMCGC collaboration; more details can be found in [^]. We fixed the 
configurations to the Coulomb gauge. 

The light quark propagators were generated using the clover fermion action at three 
values of the hopping parameter, k = 0.1385, 0.1393 and 0.1401. From an analysis of 
the light hadron spectrum j^, the second k value corresponds to a quark mass close to 
strange, where = 0.1392(1) from Mk, = 0.1394(1) from Mk- and 0.1392(1) from 
Mf, Kc = 0.1408. The 0{a) improvement term in the clover action is implemented with a 
tadpole- improved value for the clover coefficient, Csw This amounts to dividing all the gauge 
links by Uq, where we use Uq = 0.867 measured from the plaquette, and setting the coefficient 
to the tree-level value Cgw = 1- In order to improve matrix elements to the same order as 
the spectrum, we use the prescription proposed by the Fermilab group and replace the 
quark field normalisation \/2k, with ^/l — 6k, where k = uqk. 

In this simulation we truncate the NRQCD series at 0{1/Mq), where Mq is the bare 
heavy quark mass, and the action takes the form: 

S = t/j^{Dt + Ho + 6H)^ (1) 

where 

A(2) a ■ B 

Ho = and SH = -cb . (2) 

" 2Mo 2Mo ^ ^ 

Tadpole improvement of the gauge links is used throughout and the hyperfine coefficient is 
given the tree- level value cs = 1. We use the standard clover-leaf operator for the B field. 
The heavy quark propagators were computed using the evolution equation [0: 

a6H\ aHoY^ .f^ aH^Y ( ^ a5H" 



for all t, where n is the stabilising parameter. 

We generated heavy quark propagators at 10 values of (aMo,n) corresponding to (0.8,5), 
(1.0,4), (1.2,3), (1.7,2), (2.0,2), (3.0,2), (3.5,2), (4.0,2), (7.0,1) and (10.0,1). This roughly 
corresponds to a range of meson masses from Mb/3 to AMb and is sufficient for a reasonable 
investigation of heavy quark symmetry. Results in the static limit are not presented. We 
were not able to improve on the quality of the signal compared to previous results using 
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Wilson light fermions P,^, and in this case the extrapolation of the NRQCD results to the 
static limit was considered more reliable than results at the static point itself. To improve 
statistics we have also calculated the propagators on the time-reversed configurations. 

Details of the fitting analysis and extraction of the spectral quantities and decay constants 
from the heavy-light meson correlators can be found in the appendix. 



III. SYSTEMATIC ERRORS 

In this section we discuss the systematic errors associated with our choice of quark ac- 
tions and simulation parameters. In previous work we used Wilson light fermions which 
was estimated to introduce the largest systematic error into the calculation of heavy-light 
quantities. For tadpole-improved clover fermions the systematic errors are now estimated 
to be a relative error of O^^aKQcoY) ~ 4% (where we take gAqcd = ciAy = 0.185 for 
these configurations) in quantities of 0{Aqcd)- However, significant residual 0{aaAQCD) 
errors may remain in this simulation. An investigation of the scaling behaviour of quenched 
light hadron masses using a nonperturbative determination of the 0{a) improvement coeffi- 
cient [|l^] in the clover action concluded that these errors are negligible for /3"^^°^6.0 [jlT 



The remaining scaling violations at these /3's are estimated to be a few percent []TT| in the 



quenched results. In the present study tadpole- improved clover fermions are used for the 
valence quarks and staggered fermions for the dynamical quarks (which introduce O(aa^) 
errors). Since = 5.6 corresponds to roughly the same lattice spacing as that at 

pnf=o _ g_g ^ggg table m) we do not expect the associated scaling violations to be larger 
than in the quenched case. 

The truncation of the NRQCD series at 0(1/M) introduces an absolute error of 
0{AQCD{AQCD/My). This corresponds to approximately 1% errors in quantities of 
0{Aqcd) and 10% errors in the coefficients of the heavy quark expansion at Mb- In fact 
the significance of the higher order terms omitted from the action depends on the quantity 
considered. The additional terms to 0(1/M^) and the leading terms are: 

8(Mo)3 + 24Mo 16n(Mo)2 ' ^ ^ 

where V is the symmetric gauge-covariant lattice derivative and A*^^) is a discretised version 
of the continuum operator Y^Df. The first two terms, of 0(1/M^), represent the spin-orbit 
interaction (s.o.) and the Darwin term, respectively. The s.o. interaction only contributes 
to states with non-zero angular momentum and is a 1/M^ ~ 1% correction to the s.o. 
interaction arising from the light quark, which dominates (see the next section) . The Darwin 
term is a ~ 1% spin-independent shift which appears in both S and P states, but which 
does not affect spin splittings, such as the hyperfine splitting. The next order terms at 
0{1/M^) ~ 0.1% are due to corrections to the kinetic energy and discretisation corrections 
and should be even less significant. Of course, all the terms in 6Hh.o. will affect predictions 
through the change in the meson mass. However, for example the 1/M^ terms represent 
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a .1% correction to the meson mass, and hence, the shift in the meson mass has a less 
significant effect on predictions than those aheady mentioned. 

Indeed, a comparison by Ishikawa et al [|12| of the low lying S'-state heavy-light spectrum 
derived from an NRQCD action correct to 0(1/M) and that including the additional terms 
in equation ^ showed the higher order terms to be insignificant. A more detailed analysis of 
the effect of higher order terms on the S-state spectrum has been performed by Lewis and 
Woloshyn |[T^. These authors verified that the importance of the terms in the NRQCD action 
follows the expectations of naive power counting, i.e. \0{1/M^)\ < |0(1/M^)| < |0(1/M)|, 
where all 0(1/M^) terms were included. Terms of 0(1/M^) and 0(1/M^) were found to 
have a negligible effect on the S-state B meson mass splittings. 

Ishikawa et al also compared the results for the pseudoscalar and vector decay constant 
including tree-level terms to 0{1/M) and 0(1/M^) (for these quantities, corrections to 
the currents as well as to the action must be considered). The decay constants differed 
by approximately 3% at the B meson, and, 6% is, conservatively, estimated to be the 
corresponding error in fs if ^Hh.o. and the 0(1/M^) current corrections are omitted. 

An additional uncertainty is introduced through the use of the tree-level value (with 
tadpole- improvement) for the hyperfine coefficient, cb- Tadpole improvement is ex- 
pected to account for most of the renormalisation of this factor. However, an error of 
0{aKQCD{.^QCD/M)) will remain; for example this corresponds to approximately 10—30%, 
for a ~ 0.1—0.3 in the B*—B splitting, which is proportional to cb- The one-loop correction 
to Cb has been calculated [|14| and for Mq = 1.9, a bare quark mass close to Mq for this 



simulation, cb ~ 1.15—1.3, depending on the characteristic scale q* used. 

For the determination of and /s* , the factors matching the lattice currents to the full 
continuum QCD currents must be calculated. Morningstar and Shigemitsu have computed 
these factors for both the axial- vector and vector currents. The mixing at 0{a) 



among heavy-light current operators has been fully taken into account. In addition, an 
0{a{ah.QCD)) discretisation correction to the current was found and corrected for in the 
case of clover light fermions. This has been discussed in detail in 0,0, and we present an 
outline in section ^ The largest remaining perturbative error is O(a^) ~ 1—10%. 

Combining the tadpole-improved clover action for light quarks with the 0(1/M^) 
NRQCD action, Hein et al [jl8| found no significant scaling violations in results for the 
quenched B spectrum and decay constant between (3"-f^^ = 6.0 and 5.7. The results from 
a recent quenched calculation by Ishikawa et al of Jb at three (3 values in the range 
5.7—6.1 supports this finding. 

A further continuing uncertainty is the possibility of large finite volume effects. A study 
by Duncan et al in the static limit in the quenched approximation found no significant 
finite volume effects in the results for the lowest lying S-state B mesons (and in particular 
for Jb) for a box size as small as 1.3 fm. Finite volume effects are likely to be larger when 
dynamical quarks are included [^. However, since for our simulation, L = 1.6 fm, for a 
visible effect, finite volume problems would have to be very different in the dynamical case, 
despite the fairly heavy sea quark masses (m^ea is around the strange quark mass) that we 
use. A preliminary study by the SESAM collaboration [^], using sea quark masses which 
correspond to iriT^/mp = 0.69, found the finite volume effects for the p meson (which is bigger 
physically than the lowest lying S'-state B mesons) to be small, less than 5% for a lattice 
size similar to that used in this work. However, in the absence of a detailed investigation of 
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finite volume effects for B mesons at n/ = 2 this remains an unknown effect. 

Table | summarises the systematic errors expected to dominate each mass splitting and 
combination of decay constants. In general, our statistical errors are of a comparable size. 
The only significant exception is for the radial and orbital excitations of the B meson and 
the At) baryon. These particles are probably squeezed on a box of extent 1.6/m, although 
for B{1P)—B{1S) and B{2S)—B{1S) our statistical errors are also large. 

While within the B system we are close to achieving our goal of controlling the main 
systematic errors, a major uncertainty is introduced when we convert from lattice numbers to 
physical predictions. In particular, this uncertainty affects predictions through the fixing of 
the quark masses (discussed below) as well as the final conversion from lattice units to MeV. 
Ideally one would fix the lattice spacing using a quantity within the B system. However, 
with the exception of the lowest lying S'-states, the experimental results for the B spectrum 
are still preliminary and it is unfeasible to fix a, for example, from the P—S or 25" — IS* 
splitting (which have little dependence on mg). Thus, one must consider determinations of 
a from light and heavy spectroscopy; table H gives some examples for these configurations. 

As in the discussion above for heavy-light hadrons the effect of systematic errors on 
the value of the lattice spacing must be considered. We computed nip using the tadpole- 
improved clover action and results for this quantity from quenched studies are consistent 
with minimal discretisation errors in this quantity for ~ 2 GeV Similarly, the bb 



spectrum was computed with an Mv"^ NRQCD action, and in quenched studies the spin- 
averaged splittings obtained using this action show no scaling violations in the same range 
of lattice spacings 



Finite volume effects may affect rrip for a box size of (1.6 fm)^ (this is not expected to be an 
issue for heavy quarkonia). These effects are not thought to be significant for a calculation 
of rrip with similar parameters in the quenched approximation (such as that detailed in 
table |l|) p4|. However, finite volume problems are expected to be larger when dynamical 
quarks are introduced Naively, one expects mp to be overestimated, and hence 

underestimated if the meson is squeezed on a finite lattice. Although the preliminary SESAM 
study, mentioned above, suggests this may not be significant, finite volume dependence 
remains an uncertainty. 

Table || indicates there is a large discrepancy between determined from light and 
heavy spectroscopy for the HEMCGC configurations. In the quenched approximation, differ- 
ent lattice spacings are expected from quantities dominated by different physical scales and 
one must choose the most relevant quantity with which to fix a. When dynamical quarks 
are partially introduced some convergence in the estimates of is expected (although 
with rusea ~ f^strange significant residual quenching effects will probably remain) and the 
choice of quantity to determine a should be less important. However, the range in is 
comparable to that seen in the quenched approximation at (3^^^^ = 6.0, also detailed in the 
table (the same light quark and NRQCD actions were employed); IP—lS/mp = .43(1) and 
.46(1) for Uf = and 2 respectively, compared to 0.57 from experiment. One possibility is 
that introducing two flavours of sea quarks with m^ea ~ ^strange produces little effect on 
mp/T(lP— IS*). Alternatively, any effect of the sea quarks may be counteracted by more 
signiflcant flnite volume problems (or some other systematic error dependent on the number 
of sea quarks) at n/ = 2 for nip compared to = 0. 

Until this issue has been clarifled, we choose the lattice spacing from rrip to convert to 
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physical units at ri/ = 2 and use = 2.4 GeV to give a (conservative) indication of the 
uncertainty in a~^. This translates into an error in physical predictions at least as large 
as those already mentioned. We comment on how to compare our results for heavylight 
hadrons with those from quenched simulations below. 

To pin-point the B and Bs mesons from the range of meson masses simulated in this 
study we must fix the bare quark masses corresponding to the lightest quark (m, d), the 
strange quark and the bottom quark. The bare mass of the lightest quark is fixed in the 
standard way by extrapolating the pion mass to zero. As indicated above, m° is derived 
from the 0, K and K* mesons, and we use the difference in the results as an estimate of 
the uncertainty in Kg. The magnitude of this error found in the B^—B splitting and the 
ratio J'b / fss is given in table |. In order to determine M° we must first correct the meson 
simulation energies for the removal of the mass term in the NRQCD action. As detailed 
in [H we use the mass shift calculated from heavy quarkonia dispersion behaviour to obtain 
the heavy-light meson masses from the simulation energies. Thus, interpolating between our 
results for the pseudoscalar meson until we obtain the B meson mass, we find oMq = 2.1—1.8 
for = 2.0-2.4 GeV. 

In order to investigate sea quark effects we need to compare our results with similar 
calculations in the quenched approximation. While we believe the dominant systematic 
errors are under control, most of the remaining residual uncertainties are minimised if we 
compare with a simulation which uses a lattice with a similar lattice spacing and physical 
volume. From table H /3 = 6.0 in the quenched approximation roughly matches f3 = 5.6 
at n/ = 2. References and p| report on NRQCD results using tadpole-improved clover 
light quarks at this /3, on a lattice volume of ~ (1.6/m)^. A higher-order NRQCD action is 
employed for the heavy quarks (-ff^.o. is included), however, considering the negligible effect 
these terms have at the B meson mass, these quenched results are suitable for comparison. 

For the quenched simulation the lattice spacing derived from light spectroscopy is consid- 
ered most relevant for B splittings and decay constants which are expected to be dominated 
by the 'brown muck'. Thus, we initially compare the two simulations fixing the lattice spac- 
ing from rrip in both cases. The effect on the comparison of the large uncertainty in for 
n/ = 2 is also considered. 

Within our study it is premature to attempt an extrapolation of the quenched and 
unquenched data to a physically relevant number of sea quark flavours. We introduce un- 
certainties by implementing clover light valence fermions while using staggered sea fermions. 
In addition, we chirally extrapolate the valence light quarks to zero light quark mass while 
keeping the sea quark mass fixed around the strange quark mass. Thus, our simulation does 
not correspond to that for a light and heavy quark bound state with a sea of two flavours 
of physical quarks. However, inserting dynamical quarks allows the sensitivity of various 
quantities to these effects to be investigated. A more systematic approach to sea quark 
effects is reserved for later work. 

IV. SPECTRUM RESULTS 

Figure |I| presents our results for the lower lying B meson spectrum and the Af, baryon 
compared to experiment. At this initial stage our results are in broad agreement with 
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experiment. The experimental results for all but the lowest states are still uncertain. In 
particular, there are only preliminary signals for the 25" and P states. For the latter, a 
signal has been found for a Bttit resonance |^5|], which is probably has ji = 3/2, and a 
B^*^7i resonance [^, which is likely to be a superposition of various P states. Our result 
for the P state (-Bi) is obtained from an operator with quantum numbers ^Pi in the '^^'^^Lj 
nomenclature of quarkonia. Since charge conjugation is not a symmetry in the heavy-light 
system the two / = 1, J = 1 states mix and thus our operator has an overlap with both 
states. We investigated this mixing by forming a matrix of correlators with the ^Pi operator. 
However, these operators were found to be degenerate at this level of statistics. In fact, we 
are unable to resolve any splitting between the P states in this study. 

In addition, we can compare with the theoretical expectations from Heavy Quark Sym- 
metry (HQS), shown in figure ^ [^. The picture of a heavy quark surrounded by a light 
quark cloud, where the heavy quark acts merely as a colour source in the heavy quark limit, 
predicts a gross spectrum determined by the light quark degrees of freedom; this gives rise 
to meson mass splittings of 0{Aqcd) which are independent of the heavy quark at lowest 
order. In particular, the large scale features of the spectrum are due to the radial (25*) 
and orbital (IP) excitation of the light quark. At the next order are the ji = 3/2 and 1/2 
doublets of the P states, which are split due to the spin-orbit interaction. The 25'— 15 and 
P—S splittings in figure |l] are 300—600 MeV or 0{Aqcd), in agreement with this naive pic- 
ture. Similarly, heavy-light hadrons which only differ from the B meson in the light quark 
flavour, for example Bs, or the number of light quarks, for example A^, will also give rise to 
splittings independent of tjiq at lowest order. 

At finite mg, flavour and spin symmetry are broken by the kinetic energy of the heavy 
quark and the hyperfine interaction. The latter removes the degeneracy in the heavy quark 
spin, for example in the ji = 1/2 S-states and the ji = 3/2 and ji = 1/2 P-states. These 
hyperfine doublets are expected to be split by 0(Aq(^^/M) ~ 50 MeV; this compares well 
with the experimental hyperfine splitting, B*—B = 46 MeV. These splittings vanish as 1/M 
in the static limit. 

Our results are compared with experiment in more detail, taking into account the system- 
atic errors, in table |IT1|. If we assume, initially, = 2.0 GeV is a reasonable estimate of the 
lattice spacing for these configurations the only significant disagreement with experiment is 
found in the Ab—B splitting. Finite volume effects, which are not included in the estimates 
of the systematic errors, may well account for the rather high value for the A;, mass. For all 
splittings, the error from setting is at least as large as the statistical or other systematic 
errors. However, since this uncertainty leads to a positive shift in the splittings it does not 
lessen the disagreement with experiment found for the A^—B splitting. There is also less 
agreement for the Bg—B and B**—B splittings using = 2.4 GeV. 

With the mass splittings calculated for a wide range of heavy quark masses, not just 
in the region of the B meson, we are able to investigate violations of HQS at finite mq. 
The behaviour of these splittings as a function of 1/Mps is shown in figures § and ^, in 
lattice units. The E{2S)-E{1S), E{^So)^,-E{^So)^„ E{Aq)-E{1S) and E{^Pi)~E{lS) 
splittings are very weakly (and linearly) dependent on l/rriQ, consistent with behaviour 
dominated by the light quark, with small corrections due to the kinetic energy of the heavy 
quark and the hyperfine interaction; E(^So) denotes the simulation energy of a pseudoscalar 
meson, while -£'(15') denotes the spin-average of the simulation energy with the corresponding 
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vector meson. E{^Sq)k^ indicates that the hght quark mass has been interpolated to that of 
the strange quark mass. Note that by taking the spin-average the dependence of that mass 
on the hyperfine term is removed. Similarly, the hyperfine splitting in the IS" and 25* states 
depends linearly on l/rriQ over the range of masses studied, vanishing in the static limit, as 
expected theoretically. 

To quantify the heavy quark mass dependence we fit the mass splittings to the form: 

AM = Co + Ci/M + ... (5) 

where M is the pseudoscalar meson mass, and extract the coefficients Cq and Ci, which are 
expected to be of 0{Aqcd) and 0(Aq^^) respectively (with the exception of the hyperfine 
splittings where Co = 0). In fact, the coefficients can be related to the binding energy and 
expectation values of the kinetic energy of the heavy quark and the hyperfine interaction 
in the static limit for each meson or baryon. These quantities are needed in analytical 
approaches such as HQET. 

If we define the heavy-light hadron binding energy, A^*^, using: 

Mj, = mQ + A^« (6) 

where Mh is the mass of the hadron, and expand using first order perturbation theory, 

Ar =^H + ^ {H\ Ou^n \H) + ^ {H\ 0,y, \H) (7) 

where (if | and Kh represent the hadron state and the binding energy, respectively, in the 
infinite mass limit; Okin = —ip'^{D'^/2)%lj and Ohyp = — '?/'^(cr ■ B/2)ip. In principle, the 
separation of the meson mass into the heavy quark mass and binding energy in equation ^ 
must be clearly defined. However, we only consider physical mass splittings, for which the 
differences in the definitions of A^'^ cancel. This has the additional benefit that the 0(1/ M) 
coefficients of the mass splittings correspond to the physical difference in expectation values 
of the kinetic energy (and/or hyperfine interaction) of the particles in the splitting. Thus, 
we consider, for example, the 25'— IS* mass splitting of the spin-average S-states, which can 
be expressed as 

E{2S)-E{1S) = A,s - Ai5 + ^ {2S\ |25)^,^, - ^ {1S\ \1S)^,^^ (8) 

where (15*1 Okin ll'S')^^^^ denotes the physical expectation value of the kinetic energy in the B 
meson. The results are given in table ^ in physical units and compared with the theoretical 
expectation. 

Given the weak dependence of most splittings on the heavy quark mass, an estimate of 
the intercept and slope for some splittings can be made from the experimental results for the 
B and D spectrum and A;, and A^ baryons. The hyperfine splittings are strongly dependent 
on rriQ. However, the intercept is zero, and this can be used along with the splitting for the 
B meson to estimate the slope. The corresponding estimates of the intercept and slope of 
various meson and baryon mass splittings are also shown in the table. 



Table |^ shows that for quantities which are weakly dependent on l/mq, we can reliably 



extract Cq. However, it is difficult to reliably extract the slope and an increase in statistics 
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is required in order to provide quantitative predictions. Qualitatively, the results are consis- 
tent with the theoretical predictions of only small violations of HQS, 0{Aqcd/M) ~ 10%, 
around the B meson. Comparing with the estimates derived from experiment we find our 
value for the slope of the hyperfine splitting is low, reflecting the low value we obtain for 
the B*—B splitting. The absolute statistical error in the splitting is more or less constant 
with mg, while the splitting rises as the mass decreases. Thus, a low value of the slope is a 
clearer indication than the B hyperfine splitting itself, that our results are inconsistent with 
experiment. The hyperfine coefficient, c^, depends on the heavy quark mass and a determi- 
nation of this coefficient for the range of masses studied here is needed to clarify the source 
of any remaining discrepancy with experiment, for example residual quenching effects. Rea- 
sonable agreement is found between the lattice results for the E{} SQ)t,^—E{} S^jt,^ splitting 
and the estimates from experiment. This is also the case for the slope of the E{h.}j)—E{lS) 
splitting, while the intercept is too high (consistent with finite volume problems mentioned 
previously). Our results for the intercepts and slopes of spectral quantities are consistent 
with our previous estimates using Wilson light quarks detailed in reference [§]. 

A. n/ dependence 

An indication of sea quark effects can be found by comparing the Uf = 2 results with 
our previous spectrum calculation in the quenched approximation at /3 = 6.0, detailed 
in reference 0. The results for the splittings which have been calculated in both cases are 
detailed in figure |I| and table |V|, where is fixed from rrip in both cases. With the exception 
of the hyperfine splitting, the quenched simulation reproduces the experimental spectrum 
well to this level of statistical and systematic uncertainty. Furthermore, comparing with the 
Uf = 2 results, we see consistency to within l—2a. However, switching to = 2.4 GeV 
for the dynamical configurations, differences of ~ 3a appear for the Bg—B, B**—B and 
Afj—B splittings. The uncertainty in the lattice spacing must be reduced before any clear 
indications of sea quark effects can be found. 

The hyperfine splitting is a quantity where sea quark effects are expected to be seen. 
However, in practice this is not easy to investigate due to the difficulty in measuring such 
a small quantity. The quenched results for B*—B and B*—Bs, with lower statistical errors 
than at n/ = 2, are approximately half the experimental values. The splittings may be 
boosted by as much as 30% to ~ 31 and ~ 35 MeV respectively when the 1-loop corrections 
to cb are included. The statistical errors for B*—Bs are small enough that the splitting is 
still significantly below experiment, suggesting quenching effects. The corresponding nj = 2 
result is also increased to ~ 34 MeV (for = 2.0 GeV), when cb is corrected. However, 
the statistical error is large and must be reduced significantly before we can see if partially 
including dynamical quarks increases the hyperfine splitting. Using = 2.4 GeV for the 
Uf = 2 results only leads to a la difference between quenched and Uf = 2 and thus does not 
significantly change the comparison. 

V. DECAY CONSTANTS 

In the continuum, the pseudoscalar (PS) and vector (V) decay constants are defined by 
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(0| Ao \PS) = (0| g757o/i \PS) = fpsMps, (9) 
(0| Vk \Vk) = (0| q^kh \Vk) = ekfvMv, (10) 

where q and h represent 4-component light and heavy quark fields respectively. In lattice 
NRQCD, Aq and Vi are given by power series of operators in 1/M. To a given order, all 
operators with the appropriate quantum numbers and power of 1/M appear. This includes 
operators representing discretisation corrections, which vanish in the continuum. The corre- 
sponding matrix elements for these operators are combined with the renormalisation factors 
matching the lattice matrix elements to (0| Aq \PS) and (0| Vk \ V) in full QCD, where mixing 
between the operators under renormalisation must be taken into account. 

In this study we truncate the NRQCD action at 0(1/M), and match the decay constants 
to full QCD through 0{a/M). The axial- vector current is then []: 

2 

(Aq) = J2 Cjia, aM) (j^) + ©(l/M^, a^ a^, aa/M), (11) 
i=o 

where 

>f ) =(0|g757oQ|P^), (12) 



(4^)) =(0|-g757o|^g|P5), (13) 

(^P) =(0|gf^757og|P5). (14) 

q is now the light quark field in the lattice theory, and Q is related to the 2-component 
heavy quark field ip in lattice NRQCD by 

4°^ and are the tree- level operators which are obtained through the 0(1/M) inverse 
Foldy-Wouthuysen transformation 
theory. Thus, Co and Ci are 0(1), 



(2) 

Foldy-Wouthuysen transformation connecting hto Q. Ji appears at 1-loop in perturbation 



Co = 1 + apo, Ci = 1 + api, (16) 

while C2 is 0(a), 

C2 = a(p2-a), (17) 



where we use the nomenclature of reference ||T^. At this order one must also consider an 
0{aa) lattice artifact, which appears through the mixing between J^^ and jjf^ and which 



does not vanish in the infinite heavy quark mass limit |15]: 



^In the limit of zero light quark mass. 
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Jt") ={Q\aqD-ll^loQ\PS). (18) 

This factor can be viewed as a discretisation correction to Ji \ and hence, we define an 
improved operator ||T^: 

jmmp ^ j(0) ^ CAJt". (19) 

where Ca = a(l + CA/(2aMo)). Qa refiects the freedom in the definition of the improved 
operator and is cancelled by the term —aC,A appearing in C2. In practice we set it to zero, 
which one is free to do at 1-loop accuracy. 

In the lattice simulation, the number of matrix elements that must be calculated can be 
reduced by noting that 

j{disc) _ j{2) 



2Mo 

and, also, that at zero momentum on the lattice 

(4") = (^f>. (21) 

Thus, it is sufficient to compute the matrix elements corresponding to the tree-level oper- 
ators. Since these matrix elements are generated separately, the contribution to the decay 
constant from each matrix element in equation |Tl| can be analysed. For this purpose, we 
define, 

1 



f?s4^s = ^j^{Q\J?\PS). (22) 



We use fpsyMps to denote the total decay constant at both tree-level and 1-loop. 
For the vector current: 

4 

(Vk) = Cvj{a, aM) (jjf'^) + 0{1/M^, a^, a^, aa/M) (23) 



i=o 



where. 



4°^) ={0\<llkQ\Vk), (24) 

4'^) =(0|-^7.i](^^l^fc)' (25) 

= {0\ ~ q^lolkQ \V,) , (26) 

4'^) =m-^7^Q\Vk), (27) 



■2M, 







4'^) =(0|g^QI^.). (28) 
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J^^ and J^^^ correspond to tree-level operators, while the rest appear at one-loop in pertur- 
bation theory; Cyo and Cyi are 0(1), 

Cyo = 1 + aPo, Cvi = l + apX, (29) 
while Cy2 5 C*y3 and Cy4 are 0(a), 

Cv2 = a(p^ - Cv), Cv3 = aps, Cv4 = apl- (30) 

Analogous to the axial- vector case, there is a discretisation correction which can be absorbed 
into a redefinition of the zeroth order operator: 



[Jt") ={0\-aqD-l l.lkQ \Vt) . (31) 
jm^p = j(0) + CvJt", (32) 

where Cy = ct{l + Cy/(2aMo)); we set Cv = when computing Cy and Cv2- Note that, 

= -J«+2jf, (33) 
^ Jr^) = Jr, (34) 



2Mo 

and at zero momentum on the lattice, 

(4")=(4''>- (35) 

Thus, the 6 vector matrix elements can be reconstructed from the subset j = 0,1,3. For the 
purpose of analysing the individual contributions to the vector decay constant of the matrix 
elements ^4|- 12^ and we define, 

1 



fl^^^=-^{0\ji^^\V,). (36) 

We use fyy/My to denote the total decay constant at both tree-level and 1-loop. 

The matching coefficients Cj and Ca, and Cyj and Cy, have been evaluated to 1-loop 
in perturbation theory by Morningstar and Shigemitsu [p!5i|T6| . These coefficients depend 
on M and a, where we take av{q*) |0 for the strong coupling. The scale q* also depends 
on the heavy quark mass. This scale is unknown at present, although in the static limit 
Hernandez and Hill found that, q* = 2.18/a At this stage, we assume that q* for each 
mass lies somewhere between 1/a and vr/a, and compute fps and fy for these two limits. 

The coefficients Cq and C*i in the pseudoscalar case, and Cyo, Cyi, Cy2 and Cy^ in the 
vector case, include a term aln{aM) / n . When computing fs and /b* we insert aM^ for the 
heavy quark mass. There are no problems with large logarithms for the quark masses used in 
this study, however, for larger M, resumming the logarithms would be more appropriate. A 
different approach is needed when extracting the slope of the decay constants and comparing 
with the static limit. In this case, it is assumed that the meson simulated at each Mq is 
the B meson (to a better or worse approximation), and thus, oMq must be inserted in the 
logarithms for all heavy quark masses. 
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A possible concern is that the a/(aM)"' terms in the coefficients combine with nonper- 
turbative lattice errors, oc (aAgcz?)'^ in the matrix elements (^Jl^') or (•^l^'')- This, then, 
leads to contributions that look like physical {Aqcd/M)^ corrections to fs or fs*- However, 
since (Jl^) and (jk'^) are designed to match full QCD through 0{aa), these errors appear 



at 0{a'^{aAQCD)) or O^^aAQcnY)- The leading contribution to fs or fs* from these errors 
is beyond the order in the NRQCD series considered here. These issues have been discussed 
previously, in more detail, in references and [|1^ 



Theoretically, within HQET, the decay constants in the combinations, fps^/Mps and 
fv^/Mv, are non-zero and equal in the static limit, reflecting spin and flavour symmetry. 
Away from this limit, 0{1/M) corrections take the form: 



where the decay constant in the static limit, (/a/M)°^, and the coefficient of the linear 
slope, C, are expected to be of 0{Aqcd)- Thus, the ratio of the vector to pseudoscalar 
decay constant can be expressed as: 

fv4^/fps4^s = l + {Cv~ Cps)/M + .... (38) 

Another quantity to consider is the ratio of the pseudoscalar decay constants for the heavy- 
light mesons containing a u light quark and that containing a s quark. 



fpS\/ Mps\ks ( JpsV^psIks 



fpsVMpslnc \fpsV Mps\kc 



;i + (C--C-)/M + ...). (39) 



Since the heavy quark is almost a spectator in the heavy-light meson, the shift in the decay 
constant when the light quark mass is changed is expected to be roughly independent of M. 
Thus, Cps\n^ is approximately equal to Cpsln^, and the ratio should have little dependence 
on the heavy quark mass. 

The contributions to the slope of the decay constants from the 0{1/M) terms in the 
NRQCD action and corrections to the current can be identified. From ffist order perturba- 
tion theory in 1/M about the static limit: 

C = Gkin + 'idAlGhyp + dMGcorr/G, (40) 

where rfjv/ = 3 and —1 for pseudoscalar and vector mesons respectively. Thus, 

Cps — Gy = SGhyp + 2Gcorr/^- (41) 

The kinetic and hyperfine terms contribute to the decay constant through the correction to 
the meson wavef unction. 

(0|gtr(2V|P)^GH„ = 

(0| j rfyT{gtr(2)^(0), Ok^n{y)} \P)^ , (42) 
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(0| / rfyT{gtr(2)^(0), OuyM) \P)oo > (43) 
while Gcor-r is directly related to the tree-level current correction: 

(0| q^T^'H \P)^ dMG^orr/Q = (0| gtr(^)^^ \P)^ . (44) 

|P)^ represents the meson in the limit of infinite heavy quark mass and T^"^^ = 1 and a 
for the pseudoscalar and the vector, respectively. Note that equations ^2|-|^ are tree level 
expressions. This analysis is discussed in more detail in reference H. 



A. Results 

i-(0,l) 



The bare lattice matrix elements computed in the lattice simulation, fps vMps and 
fy'^'^^y/My, are given in the appendix in tables and [VI1| for the pseudoscalar, and 



tables Vlll , [TX] and ^ for the vector, for the range of masses studied. Combined with the 



appropriate matching factors, we obtain the tree-level and renormalised decay constants 
shown in tables ^ and pCll| , for the pseudoscalar and vector respectively. 

Considering, initially, oMq = 2.0, a quark mass in the region of the b quark mass (aMg = 
2.1 for = 2.0 GeV), table O shows that fpsy/Mps is reduced from the tree-level value 



by 10% for aq* = n and 16% for aq* = 1.0. Thus, the combined 1-loop corrections to the 
decay constant, 0{a), 0{aa) and 0{a/M), are significant but not larger than expected. In 
addition, the uncertainty in fs due to the unknown q* is not unduly large, ~ 6%, although 
it is bigger than the statistical errors of ~ 3%. 

We study the tree-level and renormalised currents further by analysing the individual 
contributions to the PS decay amplitude from the various operators and matching factors. 
Table pClll|, details the percentage correction of these contributions to the zeroth order tree- 



level matrix element, (^Ji^j', aMo = 2.0 and aq* = 1.0 is used for the 1-loop terms. Note 
that the contributions from Jp-* and Jf^'^ are reconstructed from (j^^'^ 



From the table we see that the 0(1/M) tree-level correction, (^J*^^-'^, reduces the decay 
constant by 13%. The 1-loop correction for this matrix element is small, ~ 1%, however, the 
other 1-loop terms, the 0{a) correction to {| J^^-*^, apo? and the 0{a/M) matrix element and 

the 0{aQ) discretisation correction, ap2 ('^'■^■') + ^ are all of a similar magnitude to 

the tree-level correction, around 10%. However, none of the terms are unduly large and there 
is no indication that NRQCD is breaking down. Since the terms contribute with different 
signs there is some cancellation between them and the combined 1-loop correction is —14% 
for aq* = 1.0, the same size as the 0(1/M) tree-level correction. Thus, the total correction 
to /?]v/M^ is -27%; this falls to -21% for aq* = vr, of which -8% is due to the 1-loop 
terms. 



The picture is fairly similar for the vector decay constant. Table |X11| shows fy^My at 



aMo = 2.0 is reduced from the tree-level value by 14% and 22% when the matching factors 
are included with ag* = vr and 1.0, respectively; the dependence on (f is moderate. In the 
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tree-level case, fpsV^PS is slightly below fvy/My, while for the same q* they are roughly 
equal. However, there is no reason for the characteristic scale for the vector to be the same 
as that for the pseudoscalar, and thus, it is difficult to compare the two at this stage. 

The contributions to the vector decay constant from the individual matrix elements and 
matching coefficients are also detailed in table pClll| . We see that the 1-loop terms, apo and 
ap2 ^J*-^-*^ + a l^J'^^^^'^ are all the same order of magnitude as in the PS case, and combine to 
give a —20% correction to the zeroth order vector current. The remaining 1-loop terms, all 
0{a/M), are very small, of order 1%. However, since the 0(1/M) tree-level term is smaller 
and with opposite sign, compared to the PS case, the overall correction to ('^'■°'') is smaller: 
— 17% (—10%) for aq* = 1.0 (tt). The relative size of the individual corrections compared 
to those for the PS is discussed later. 

Interpolating the results for fpsV^PS and fyy/My to Mq = Mq, and converting into 
physical units we obtain the values for Jb and fp* shown in table PC1V| ; the results for fs^ and 



fsJfB are also given. Table pClV| presents the predictions for the decay constants derived 
using both = 2.0 and 2.4 GeV. We take the difference in the decay constants for the 
different choices of lattice spacing as an estimate of the error in due to uncertainty in 
the scale. 

Combined with estimates of the other dominant systematic errors, detailed in table |, 
our predictions for the pseudoscalar and vector decay constants are: 

Jb = 186(5) (19) (9) (13) (+50) MeV (45) 

fB, = 215(3)(22)(9)(15)(+49)(-5) MeV (46) 

fB* = 181(6)(18)(9)(13)(+55) MeV (47) 

fB* = 213(4)(21)(9)(15)(+60)(-4) MeV (48) 

The first error is statistical. The second is the estimate of the O(a^) perturbative error 
~ 10%. The third and fifth errors are the estimates of uncertainty arising from the light 
quark discretisation errors and the uncertainty in a~^. Note that the latter includes the 
corresponding change in aMg. For Jb^ there is an additional error due to the uncertainty in 

Kg- 

The fourth error is due to the truncation of the NRQCD series. To estimate this, one 
must consider the next order corrections, 0(1/M^), to the current and to the NRQCD 
action (given in equation p. Naively, one expects these terms to be of order ~ 1%. As 
mentioned in section |III| Ishikawa et al [|12| performed a tree-level analysis of the 0(1/M^) 



terms for fps and /y. These authors found that in the region of the B meson SHh.o. had no 
effect on the PS or V matrix elements. For lighter values of qMq the higher order terms in 
the action tended to increase the matrix elements. 
The individual tree-level current corrections are, 

= ^^-rz.^^' (49) 
J^'^ = ■ BQ, (50) 

J^'^ = -^qT2iga-EQ, (51) 
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— * 

where a = 707 and S = diag{B,a). T = 7570 and 7/,. for the pseudoscalar and vector 
respectively. These terms were each found to be less than 2% in magnitude. The corrections 
contribute with differing signs and there is a cancellation between them which leads to a 
small overall decrease in /ba/Mb of ~ 3%. Since at 0(1/M) the 1-loop terms are as large as 
those at tree-level, as discussed previously, we estimate the overall uncertainty in the decay 
constants due to the omitted 0(1/M^) and 0(a/M^) terms as 6%. 

In total, adding the statistical and systematic errors in quadrature, there is a 14% un- 
certainty in the decay constants. This is much smaller than the ~ 30% error due to the 
uncertainty in a~^. Thus, a determination of of around ~ 10% is possible, but only if 
the systematic error in is reduced significantly. 

In contrast, fssl fs and fs'/fB* are quantities for which most systematic errors cancel 
out. These ratios are independent of the heavy quark mass (see below), while table |XIV| 
indicates that it is insensitive to the 1-loop corrections. Furthermore, fsj fs and fB*/ fs* 
are not dependent on the lattice spacing. The only significant error is that due to the 
uncertainty in Kg- We find, 

fsJfB = 1.14(2)(-2), (52) 
fs^/fB* = 1.14(2)(-2) (53) 

where the first error is statistical and the second is due to k^. 

We now study the dependence of fpsV^PS on the light and heavy quark mass. Con- 
sidering the light quark mass dependence initially, we find the decay constant is linearly 
dependent on rriq as shown in an example of the chiral extrapolation of this quantity in the 
appendix for oMq = 1.0. As discussed in section |V| the dependence on anig is expected to 
be roughly the same for all aMg, i.e. Cps is insensitive to the light quark mass. Figure |^ 
gives the tree-level values for fpsV^PS for all heavy quark masses and ki. 

We perform a fit to the decay constant for fixed k/ using a functional form motivated by 
equation P7| , 

fps^/M^s = Co + Ci/M + ... (54) 

where Co = {fVM)°^ and Ci/Cq = CpsIki- The slopes for each ki and the chiral limit 
are given in table |XV| . We find the slope for the lighter two k's (corresponding to rriq ~ 
rus and smaller) are consistent with the chiral limit, Cp5 ~ 1 in lattice units. However, 
from the results for ki = 0.1385, we see there is a much stronger dependence of the decay 
constant on l/{aMps) for ruq heavier than m^, although with such large errors it is only 
a 2(7 effect. This behaviour may be due to discretisation errors in Cp5 (which is a light 
quark quantity) arising from the light quark action, which increase with the mass of the 
light quark: 0{aamg) and 0((amg)^) for the clover action. Since Cpsl^s ~ C'psUc) f^e ratio 
{fpsy/Mps)i^J{fps\/Mps)K,, has almost no dependence on the heavy quark mass, as seen 
in figure ^ (b). 

The heavy quark mass dependence of fpsVMps at tree-level and 1-loop is compared 
in figure ^ for ki = Kc- In the tree-level case there is a steep dependence on the heavy 
quark mass, and this, initially suggested higher orders in the heavy quark expansion would 
be needed at the B meson 0; fpsV^PS is reduced from the static limit by ~ 35% at 
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the B meson. However, figure ^ shows the slope is dramatically reduced when the 1-loop 
corrections are included; fpsV^PS at the B meson is reduced by ~ 17% compared to the 
static limit for aq* = tt, while there is almost no dependence on M for aq* = 1.0. 



We perform a fit to the data sets using equation the results are detailed in table pCVl 



The relative slope at tree-level is ~ —1 in lattice units or ~ —2 GeV, using a^^ = 2.0 GeV. 
This is well above the naive expectation of —300 to —500 MeV {0{Aqcd))- However, at 
1-loop this is reduced to anywhere between —1 GeV and 0, depending on q*. Note that 
q* depends on the mass, and this dependence must be determined before the slope can 
be extracted correctly. However, the results with fixed q* provide a rough bound on the 
size of the slope. Previous calculations in the quenched approximation using the clover 
action and the Fermilab approach find the slope to be close to the upper limit around 



— 1 GeV p0|-|33|. Note that the values of q* are also not known for this method. 

There are clear advantages to simulating at finite heavy quark mass with an inverse 
lattice spacing around 2.0 GeV. The combined 1-loop corrections grow as M tends towards 
the static limit (with the a In(aM) / vr terms in the matching coefficients fixed to a In(aMQ)), 
as does the accompanying uncertainty due to the omission of higher orders in the matching 
coefficients. Thus, for a coarser lattice spacing, i.e. with a larger value for oMq, (and 
indeed for the static limit) a nonperturbative determination of the matching factors becomes 
essential. Certainly this is required for a determination of the slope of the decay constant. 
For a finer lattice, oMq becomes smaller, and at some point NRQCD breaks down: the bare 
lattice matrix elements grow larger, i.e. increasingly higher orders in the NRQCD expansion 
become important, while the a/{aM) terms in the matching coefficients also become large 
and perturbation theory is not well controlled. 

The interplay between the individual PS matrix elements and their effect on the slope of 
the decay constant can be seen in figure |^. The tree-level and 1-loop corrections to the zeroth 
order matrix element are shown (as percentages) as a function of l/(aMo) for aq* = 1.0. 
The 0(1/M) and 0{a/M) terms, ('^^^■') and ap2 ('^^^■') and api (^J'^^^'^, respectively, vanish 
in the static limit and grow as uMq becomes small. The discretisation correction has very 
little dependence on the heavy quark mass, while the 0{a) correction to (<^''°'*) decreases 
dramatically from the static limit, and goes through zero as oMq becomes lighter. 

The corrections contribute with differing signs leading to some cancellation between 
them. Towards the heavy mass limit, l/(aMo) < 0.5, the 0{1/M) and 0{a/M) terms 
approximately cancel each other. This leaves the 0{a) and 0{aa) terms which are both 
large in this region, ~ —20% each, leading to a very large 1-loop correction to fpgy/Mps- 
In contrast, for the lighter values of qMq < 1.2, apo is much smaller, while the 0{a/M) 
terms are now large enough to approximately cancel the discretisation correction. Thus, 
surprisingly, at 1/aMps ~ 0.6, the decay constant is roughly equal to the tree-level value. 
Note, however, that in this region the 1-loop terms ap2 (^''^^) and (^J'^*'^'^^ are not small at 
around 10%. 

Thus, the almost complete removal of any heavy quark mass dependence when the 1-loop 
terms are included using aq* = 1.0, is due to the large 1-loop correction in the large mass 
limit (where 0(1/M) terms are small), ~ —40% in the static limit, and a large tree-level 
0(1/M) correction for the lighter meson masses (where the 1-loop terms cancel), ~ —30% 
for aMo = 0.8. Since the tree-level correction is so large for the lighter meson masses, higher 
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orders in the NRQCD expansion should be considered in this region. 

The vector decay constant behaves in a similar way as a function of the heavy quark 
mass to the pseudoscalar. Figure || details the individual percentage corrections to fy^ ^ My 



versus l/(aMo). The bare matrix elements, (^Jl y , i = 1,4 are related to each other as given 
in equations |33l- 13^. In the heavy quark limit, these expressions simplify and furthermore 
one finds additional relations between the vector and pseudoscalar current corrections. If 
we re-write jj^^ and jj^^ in the form, 

4'^ = (-7± ■ + ikQ (55) 

4^^ = -^j^^qikDk-fkQ (56) 

and take \ Vk) = [Q^jk^] |0), then by assuming that the contraction of QQ^ is spin diagonal (in 
the large M limit) we can deduce, 

{0\qi-fiD,)^kQ\Vk) it^k) = {0\qi^kDk)-fkQ\Vk), (57) 

i.e. (4^"*) — ~ (4^'*)- Tables |I^ and ^ show that our results agree with this relation for 
the heavier values of M>3.0. Thus, equation ^ reduces to (4^'') = 3 (4^^) = ~3 (4^'')- 
Following the same argument, and taking IPS') = [(5^75^] |0), one can show that (^J^^ 
for the PS is equal to (4^'') for the V, 

(0| qhj D,)lolkQ \Vk) = (0| g(7, B,)757oQ \PS) (58) 

for any j. Thus the corresponding discretisation corrections are also equal, i^J^'^^'^'^^ = 
j{disc)\ (Jl^') = (Jl^') = —^(jk^^y, the latter clearly comes from spin-symmetry. 



The results in figures and ^ in the large mass limit agree with these expressions. In 



general the terms which spoil equations q3 aiid|5g at finite M are due to the 0{1/M) terms 



in the action modifying the meson wavefunction from the static limit (as in equation ^ 



and ^for the PS decay constant), and are small since they appear as an 0(1/M) correction 



to the matrix element. 



B. Making contact with HQET 

As discussed in section 0, one can attempt to analyse the 1 /M corrections to the static 
limit within the framework of HQET. Equations ^ and suggest that by taking appropriate 
combinations of fpsy/Mps and fyy/My, Gkin and Ghyp can be separated; since (^J^^^^ is 
calculated separately Gcorr can easily be obtained. In particular, the spin average of the PS 
and V decay constants (without <^J*^^^^), /v^M, cancels the hyperfine contribution, and the 
slope of this quantity is purely determined by Gkin- Conversely, the ratio of these quantities 
cancels Gkin and a determination of the slope gives Ghyp- 

Note that the renormalisation factors calculated by Morningstar and Shigemitsu are not 
those required for a detailed investigation of spin and flavour symmetry within HQET. In 
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references [|T5[ and the lattice NRQCD matrix elements are matched directly to full 
QCD while, the factors matching lattice NRQCD to continuum HQET are the ones that 
are needed. In particular, only the tree-level decay constants will obey the spin symmetry 
relation that the ratio of V to PS is 1 in the static limit. Thus, we restrict ourselves to a 
tree-level analysis of spin symmetry and the origin of the slope C. One-loop matching to 
full QCD would introduce a short distance correction factor to the ratio fvV M^v / fpsV^PS 
of Vstatic = [1 - tt(8/(37r))]. 



Figure ^ presents the results for fvM, compared to the PS decay constant, with and 
without the tree- level 0(1/M) current correction. In the heavy mass limit, all three quanti- 
ties converge indicating the dominant contribution to the tree- level Cps is from Gkin- The 
figure also shows fpsVMps/ fv\/Mv- Without ('^^^■') this ratio is positive (indicating Ghyp 
is positive) and tends to 1 in the static limit, consistent with spin-symmetry; when the cur- 
rent correction is included fvV^v < fpsV^PS- Gcorr can be very accurately determined 
by extrapolating 2Mq (^J'^y to the static limit, displayed in figure |^ (c). 

We perform fits to these combinations of PS and V decay constants and extract Gkin, 
Ghyp and Gcorr- The results are given in table pCVll] . Note that only the slopes of physical 
combinations of the decay constants, i.e. fps^/Mps, fy/M and fpsy/Mps/ fyy/My, are 
expected to be 0{aKQCD)- The ratio of decay constants is roughly in agreement with this 
picture; the slope of the ratio is small, while Gcorr and Ghyp are individually much larger. 
However, the physical decay constants themselves, dominated by Gkin-, have a much larger 
slope. Considering the strong dependence on M of the 0{a) correction to fpl^/Mps seen 
in the results of Morningstar and Shigemitsu, Gkm is likely to be significantly affected by 
renormalisation. Similarly, the 1-loop corrections to Gcorr are likely to be large. 

The results in table PCVII| through the use of the clover action for the light quarks and 



extrapolating to the chiral limit, are an improvement on our previous results with Wilson 
light fermions published in reference 0. Numerically, the difference between the results 
is mainly due to the extrapolation to in the clover case since at tree-level the decay 
constant does not change significantly between the two light quark actions. 

C. Uf dependence 

We compare our results with those calculated as part of the 0(1/M^) quenched simula- 
tion at = 6.0 already mentioned in section |IV A| . In reference we quote, 

/B = 147(ll)(^L)(9)(6)MeV (59) 
fs. = 175(08)(lIo)(ll)(7)(io') MeV (60) 
fBjfB = l.20mtt) (61) 

for the quenched result. The errors are calculated in a similar way to those in this paper and, 
with the exception of the uncertainty due to fixing a~^, are of a similar size. The first error 



^The Wilson results in reference Q correspond to a light quark mass greater than and thus 
the slope of the decay constant is much larger than in the chiral limit, as discussed previously. 
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corresponds to statistical errors and those errors due to extrapolation/interpolation in k and 
oMq. The second error indicates the uncertainty in calculated using = 2.0 GeV and 
1.8 GeV. This range represents the spread in derived from various light spectroscopic 
quantities. Higher order perturbative and relativistic uncertainties are given by the third 
error and the fourth error is due to discretisation corrections. For fs, the error in fixing Kg 
is found by fixing the strange quark mass from the K meson (for the central value) and the 
meson; the K* meson gave the same value of Hs as the (p meson. 

The quenched calculation includes 6Hh o. in the action and the tree-level 0(1/M^) current 
corrections. Since the 0(1/M^) current corrections were calculated separately we can easily 
omit them in order to make a better comparison with our present 0{1/M) nf = 2 results. 
Apart from this the quenched simulation and analysis remains the same as detailed in 
reference 0. The modified quenched results are presented in table PCV111| . As mentioned in 



section |V A| , the contributions from SHh o have a minimal effect on the pseudoscalar matrix 
elements in the region of the B meson. Thus, there should only be a very small uncertainty 
in the comparison between quenched and nj = 2 results due to the differing actions. In fact, 
since SHh.o. tends to increase Jb, removing it would increase rather than decrease the sea 
quark effect that we see. 

The comparison is made at tree- level, and at 1-loop for aq* = tt and aq* = 1.0 in 
table PCV111| . In perturbation theory the effects of sea quarks appear beyond 1-loop. Thus, 



aq* is likely to be very similar for the same Ma values for quenched and Uf = 2. Since the 
values of a and hence oMq are closely matched for the two simulations, it is reasonable to 
compare at fixed aq*. 

Considering, in addition, that the quenched and unquenched simulations are very similar 
in method and analysis we expect the systematic errors to be similar (and correlated) in 
both cases. Thus, initially, we make a comparison considering only the statistical errors. 
Table PCVIII| shows that fs is a quantity which is sensitive to the presence of sea quarks 



and that the quenched value for fs may be substantially lower than the value in full QCD. 
We find a 23—25% or > Sa increase in Jb, when two flavours of sea quarks are included. 
The increase is roughly the same at tree-level and at 1-loop for aq* = vr and aq* = 1.0. If 
the actual value of aq* for Uf = turns out to be higher than that for nj = 2 then the sea 
quark effect becomes smaller. In the reverse case the effect would be larger. 

Alternatively, we can ignore the fact that the systematic errors between the quenched and 
unquenched results are probably highly correlated. Combining the statistical and systematic 
errors in quadrature, we obtain fs = 149(19) MeV at Uf = 0, compared to 186(25) MeV 
at Uf = 2; the lattice spacing is fixed using rup in both cases, i.e. the uncertainty arising 
from fixing a is ignored, and the central value is obtained from the average of the results for 
aq* = 1.0 and vr. In this case, the 27 MeV shift in fs is a la effect. Note that, considering 
the uncertainty in a~^ for Uf = 2, using a~^ = 2.4 GeV will lead to a larger sea quark effect. 
This increase includes any effect due to finite volume problems for the determination of a~^ 



from nip at nj = 2. As mentioned in section |T| a is expected to be underestimated if the 



^The significaiice of the shift in the decay constant is found by comparing the difference between 
/b at = 2 and with the statistical errors added in quadrature. 
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p meson is squeezed. 

Comparing the central values for /^^ the increase with n/ is slightly smaller, ~ 20%, but 
statistically more significant at between 5 and 6.5(T. The error in /^^ due to the uncertainty 
in Ks clouds the comparison for this quantity slightly. For the Uf = 2 configurations, 
and rriK give consistent values for Kg, while ttlk* gives a slightly higher value (see section y). 
In contrast, for the quenched results the error in Kg is larger, and as mentioned previously, 
is computed using from using ttik and m^; ttlk* and give consistent results. Thus, if 

is used the quenched decay constant is increased by 7 MeV, while the Uf = 2 value is 
unchanged. The difference between the two values for fs^ is reduced to 14—16% or a 4— 5.5a 
effect (a similar decrease is found using mx*)- A detailed study of how to fix Kg is needed 
to improve the comparison of quenched and unquenched results for Jb^- If the systematic 
errors are taken in account, fs^ = 180(20) MeV at = and 215(29) MeV for n/ = 2; the 
central value is found using Kg from m^, but the associated uncertainty in Kg is included in 
the error. Assuming the errors are independent, the 35 MeV shift corresponds to la. 

The effect of sea quarks in the decay constants have also been studied by the MILC 
collaboration [Q. fs was found to increase by ~ 13% from the quenched to the Uf = 2 
result. The Wilson action was employed for the heavy and light quark at several (3 values 
at both nf = and 2. Extrapolating the results to the continuum limit these authors find 
/s = 159(11) (stat) (+22 - 9){sys){+21 - 0){quench) MeV. Clearly systematic errors (of 
which discretisation errors are the dominant source) are as large an effect as quenching. In 
this situation the size of the sea quark effects is very sensitive to the interpolation of the 
Hf = results to a lattice spacing which coincides with that found at = 2 to enable a 
comparison. However, until the systematic errors are explored thoroughly at = 2, the 
difference between our results and those from the MILC collaboration are not significant. 

Since the matching factors have been calculated over a wide range of M, we can also 
investigate the mass dependence of the sea quark effects. Figure presents the results for 
tree-level and 1-loop for two values of aq*. The roughly constant shift between nf = and 
Uf = 2 for the whole range of meson masses from ~ Mb /3 to AMb indicates the sea quark 
effects are associated with the light quark rather than dependent on the heavy quark. We 
also compare the results for the ratio /bs/ fs- The large sea quark effects we see in /^^ and 
Jb individually appear to cancel in the ratio. There is only a ~ 5% or 1.3cr decrease in 
IbsI Ib when two flavours of sea quarks are introduced. If Kg from is used the results 
are consistent to within less than 0.5(T. 



VI. CONCLUSIONS 

We presented comprehensive results for the spectrum and decay constants of the B meson 
including the effects of two flavours of dynamical quarks. In addition, an investigation of the 
effects of sea quarks through a comparison with quenched simulations was performed. At 
present, we do not attempt to extrapolate in but look for initial indications of quenching 
effects. 

Results, summarised in table |ITT| , were presented for the lower lying 5* and P meson 
states and the A^ baryon. With the exception of the Ab—B splitting, the estimates of 
the systematic errors for all mass splittings were comparable to or less than the statistical 
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errors. The Ab—B splitting represented the only significant disagreement with experiment. 
This discrepancy is probably due to finite volume effects. 

A detailed analysis of the heavy quark mass dependence of physical spectral mass split- 
tings was performed. We found the behaviour with M to fulfill the naive expectation: 
0{Aqcd) quantities are weakly dependent on M, while spin-dependent splittings vanish as 
1/M in the large mass limit. Estimates of the slope and intercept of these quantities were 
extracted and agreement was found with the theoretical expectation from HQS. A compar- 
ison was also made with experiment. This showed that the slope for the hyperfine splitting 
was too low in our simulation. Residual quenching effects are a possible explanation for this 
discrepancy. However, we have underestimated the splitting through the use of the tree-level 
value of the coefficient of the hyperfine term in our simulation (this underestimate is not 
significant within statistics at the B meson), cb depends on M and improved estimates of 
this coefficient must be included for all qMq before the origin of the discrepancy is clear. 

We presented a comparison of our results with a similar quenched calculation at (3^^^^ = 
6.0 using a 0(1/M^) action detailed in [0. A comparison with a higher order NRQCD 
action is possible as the additional terms in the action do not significantly affect the mass 
splittings considered. The quenched results reproduce the experimental spectrum well to 
the present level of accuracy, and there is no significant difference between the quenched 
and Hf = 2 results. 

In addition, results were presented for the /b, /b^, /b* and fs* decay constants at tree- 
level and fully consistent to 1-loop in perturbation theory. We found the combined 1-loop 
corrections are moderate and reduce the tree-level value of Jb by 10—16% depending on q*. 
Our predictions for the PS and V decay constants are 

fB = 186(5) (stat) (19) (pert) (9) (disc) (13) (NRQCD) (+50) (a-^) MeV, 
fB, = 215(3) (stat) (22) (pert) (9) (disc) (15) (NRQCD) (+49) (a^^) (-5) (/€,) MeV, 

= 181 (6) (stat) (18) (pert) (9) (disc) (13) (NRQCD) (+55) (a-i) MeV, 
/b. = 213(4) (stat) (21) (pert) (9) (disc) (15) (NRQCD) (+60) (a~^) (-4) (K,)MeV. 

Adding the statistical and systematic errors in quadrature (errors 1—4), there is a combined 
uncertainty of 14% in the decay constants. This indicates a determination of Jb to around 
10% is possible using an 0{1/M) NRQCD action with clover light fermions. However, first 
of all the error due to the uncertainty in the lattice spacing (~ 30%) must be significantly 
reduced. Conversely, for fss/ Ib and fsi/fB* most systematic errors cancel, and these 
quantities can be more easily determined accurately. We find 

/B7/B = l-14(2)(stat)(-2)(/c,), 
=1.14(2)(stat)(-2)(/€,). 

We presented a detailed study of the heavy quark mass dependence of the decay con- 
stants. We find the steep dependence of the tree-level PS decay constant, for which 
Cps ~ 2 GeV, is dramatically reduced when the 1-loop terms are included; Cps = 1—0 GeV 
depending on q* . We show how this occurs through the cancellation between the various 
1-loop and tree-level current corrections as the individual terms vary with the heavy quark 
mass. 
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We have shown that simulating directly at the B with such that oMq is in the region 
of 2 is advantageous. In particular, this value of qMq sits in the window between large 
perturbative corrections for larger oMq and significant contributions from higher orders in 
NRQCD for smaller aMg. Nevertheless, the 1-loop matching factors should be checked by a 
nonperturbative determination; this is essential for a determination of the slope. Similarly, 
for the smaller meson masses studied here, higher orders in the NRQCD expansion should 
be considered. However, there are no indications that NRQCD is breaking down for the 
range of masses studied. 

In addition, we present a tree- level analysis of the origin of the slope Cps, and separate 
the contributions from the hyperfine and kinetic energy terms in the NRQCD action and 
the tree-level current correction. We find Gkin ~ 1-8 GeV dominates Cps, while the slope 
of the ratio fpsV^^Ps/ fvV^v, related to Ghyp and Gcorr is 0{Aqcd), consistent with 
naive expectation from HQS. However, considering the results for the 1-loop corrections 
discussed previously (which are not suitable for a study of HQS), Gkin is likely to be reduced 
significantly by renormalisation and 1-loop corrections to Gcorr are also likely to be large. 

With this extensive analysis of the decay constants we are able to perform a detailed 
comparison of our results with those from the 0(1/M^) quenched simulation at /5 = 6.0. 
We find that fs appears to be considerably larger at n/ = 2 compared to that at n/ = 0. 
Clearly, fs is a quantity which is sensitive to internal quark loops, and further work is 
needed to investigate the dependence of the decay constant on the mass and the number of 
flavours of sea quarks. We find the sensitivity to nj cancels in the ratio fss/ fs- 
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VIII. APPENDIX 

In this section we illustrate the fitting analysis and extraction of the masses and am- 
plitudes of the heavy-light meson correlators. The fitting method is described in detail in 
reference 0, and since the data is of similar quality in terms of efficacy of the smearing 
functions etc. we only present additional fits not shown in [§] and [§. 

A. Spectrum 

Considering the spectrum initially, the simplest quantities to extract are the ^Sq and 
^5*1 ground state energies, E(}So) and E{^Si) respectively. A multiple exponential 'vector' 
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fit was performed to tlie C(/, 1) and C{1,2) correlators, wliere C(/,l) (C(/,2)) denotes a 
correlator with the heavy quark smeared at the source with a hydrogen-like ground (first 
excited) state wavefunction and local at the sink. The operators inserted at the source and 



sink constructed from heavy and light quarks are given in table PCIX . 

An example of the quality and stability of the fits is shown in figure |Tl| for aMo = 1.0 
and Hi = 0.1385. The ground state energy extracted from a single exponential fit [uexp = 1) 
in the range tmin > 12 is in agreement with E{^So) from a Uexp = 2 fit for 3 < < 10, 
and we feel confident that we have minimal excited state contamination. The excited state 
energy is stable with tmin, however, a rig^p = 3 fit is needed to confirm this state. Similar 
results are found for all Mo and ki and we chose a fitting range of 5 — 20 with rie^p = 2 for 
both ^5*0 and ^Si mesons. 

Tables pCX| and pCXll| detail the energies extracted. The results are chirally extrapolated 
to zero light quark mass and we find only linear dependence in 1/ki; figure |12| shows the 
chiral extrapolation of E^-^ for oMq = 1.0. 

The hyperfine splittings, E(^Si)—E(^So) and E{2^Si)—E{2^So), are calculated from the 
differences of 100 bootstrap samples of the Uexp = 2 fits. In order to study the dependence of 
these and other splittings on the heavy quark mass, the meson mass [Mps) is calculated using 
the mass shifts, A = Mps — E{^So), obtained from heavy quarkonia dispersion behaviour, 
given in reference The mass shifts are given in table pCX[ 

The singlet and triplet / = 1 states, ^Pi, ^Pi, ^P2 and 3Po, were also computed for 
aMo = 1.0, 2.0 and 4.0. The corresponding operators used in the simulation are given in 



table PCIX| . As noted in we cannot resolve the spin splittings between these states with 



the statistics available. Thus, we present the results for the ^Pi state and ignore the mixing 
with the ^Pi state. Figure 13 shows a riexp = 1 fit to the C{1, 2) correlator is stable as 



is varied and we chose the fitting range 4 — 10 for all aMg and k;. The ground and first 
excited state smearing functions did provide significantly different overlap with the lowest 
two states. However, since the signal died out around t ~ 10 there were too many higher 



states contributing to enable a 1, 2 or 3 exponential fit to C(/, 1) and C{1,2). Table pCXIV 
gives the results for different oMq and ki. 

In addition to the meson spectrum the = baryon containing one heavy quark, Aq, 
was computed. The lattice operator, 

ea6cEQa(^i)E(?')'^(^2)C757og^(x2)0(|xi (62) 

Xl X2 

has the correct quantum numbers for this particle, where a, b and c are the colour indices 
and a is the spinor index. represents the (hydrogen-like) smearing function and C is the 
charge conjugation matrix 7072- Figure presents the energy of the ground state extracted 
from multi-exponential fits to C{1, 1) for qMq and ki = 0.1385. E{Aq) is stable with t^m 
for Hexp = 1- This is also the case when including another exponential in fitting to the same 
correlator. We found the fitting range 7 — 20 for the n = 1 fit is sufficient for all aMo and 



Ki, and the corresponding values of the ground state energy are given in table PCXV 
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B. Decay Constants 



The zeroth order pseudoscalar and vector decay constants, fp^y/Mps and /y ^ a/ My 
respectively, were extracted in the standard way from simultaneously fitting to the C{1, 1) 
and C(l, 1) correlators. The amplitudes of these correlators are related to the decay constant 
via: 

C{1, 1) = ZiZ,e'^\ (63) 
C(l, 1) = Z^e"^*. (64) 

where, 

V2Zi = /(°)Vm. (65) 
Note that the correction to the current is not included in C(/, 1). The typical quality of the 



fits is illustrated in figure ^ for the PS meson. We found 7 — 20 to be the optimal fitting 
range for PS and V mesons for all oMq and ki. The corresponding amplitudes obtained 
from 100 bootstrap samples of the fit parameters are presented in tables and |VII1] for the 
PS and V respectively. The results are extrapolated to zero light quark mass, an example 



of which is shown in figure for the PS decay constant. 

The tree-level corrections to the PS and V currents, fpgy/Mps and fy^y/My respec- 
tively, are obtained separately by computing the jackknife ratio of the correlator with the 
tree-level current correction operator inserted at the sink (C(/j(i),l)) with 1); J*^^^ is 
given in equations |1^ and ^ for the PS and V respectively. In the limit of large times this 
ratio tends to f^^^ ^/M / f^^y/Jd . From fi gure ^ we see that 8 — 20 is a reasonable fitting 



range, and this was found to be optimal for all oMq and ki. Tables |V11| and ^ give the 



resulting tree-level corrections to the decay constant for the PS and V mesons respectively. 
The only 1-loop correction which must be calculated is fy^y/My- This is obtained in the 
same way as the tree-level correction. The results are detailed in table ^ 
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TABLES 



TABLE L The systematic errors expected to dominate the predictions of mass sphttings and 
decay constants of the B meson. 



Quantity 


Dominant Errors 


Size 


B*-B 


Pert, correction to 


0{as{kQCDlM)) ~ 10-30% 


Bs-B 


Determination of Kg 
Light quark discretisation 
Finite volume 


~ 10% 

0(iakQCD?), 0(a{akQCD)) < 5% 
5% ? 


B{2S)-B{1S) 
B** -B 
Ab-B 


Finite volume 
Light quark discretisation 


< 10% ? 

OiiaAQCDf), OiaiaAQcn)) < 5% 


Jb, Ib* 


Light quark discretisation 
Pert, matching factors 
Truncation of NRQCD series 
Finite volume 


OiiaAgcD?) < 5% 
0(a2) ~ 1-10% 
0(1/M2) ~ 6% 
5% ? 


!bI!b, 


Determination of Kg 


-2% 


TABLE IL Values of the inverse lattice spacing in GeV as determined from light hadrons and 
heavy quarkonia [23 for this ensemble and a similar quenched ensemble at f5^i^^ = 6.0 [23, |]. 




P 


r{lP-lS) T{2S-1S) 


/3"/=2 = 5.6 
/?"/=o = 6.0 


1.97(3) 
1.93(3) 


2.44(7) 2.37(10) 
2.59(5) 2.45(8) 
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TABLE III. Predictions for various mass splittings compared to experiment. The results are 
converted into physical units (MeV) using = 2.0 GeV. The first error shown is statistical while 
the second is the estimate of the systematic errors which we expect to be dominant, as detailed 
in table |[ the third indicates the change in the prediction if a^^ = 2.4 GeV is used and includes 
the change due to the change in Mq. For Bs—B the systematic errors are split up between the 
light quark discretisation error and that due to the uncertainty in Kg- The systematic error for the 
hyperfine splitting indicates the change in the central value if is increased by 30%. 





/3"/=2 = 5.6 


expt 


B* - B 


20(12)(+6)(+6) 


45.7(4) 


Bt - Bs 


26(8)(+10)(+9) 


47(3) 


Bs-B 


98(8)(5)(-12)(+20) 


90(2) 


B{2S) - B{IS) 


540(120)(23)(+110) 


~ 580 


B,{2S) - B,{1S) 


500(80)(20)(+100) 




B*{2S) - B{2S) 


78(32)(10)(+22) 




B** - B 


520(40)(20)(+110) 


350-500 


Ab-B 


560(40)(22)(+110) 


345(9) 
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TABLE IV. The coefficients of various spectral quantities, as determined from first order per- 
turbation tlieory in 1/M. The results are in physical units, GeV and GeV^ for Co and Ci respec- 
tively, converted using = 2.0 GeV. The hyperfine term has been ignored for the ^Pi state. 
Unless stated otherwise, ni = Kc- 



Quantity 



Coefficient 



Results 



Expectation 



Expt 



E{2S) 


-E{IS) 


Co = 


A25 - 


-Ms 


0.38(10) 


+Aqcd 








Ci — (25 Okin 








; 








\'^^)phys 


0.60(40) 


H-ve 










Co 


= 


0.002(8) 





~ 






Ci = 


{is\o,y,\is);f 


0.08(4) 




~ 0.28 




-E{2^So).^ 




Co 


= 


-0.02(2) 





- 






Ci = 


(25| Ohyp\2S)p'f^ 


0.56(28) 




- 




- Ei^So)^, 


c 


- 






0.106(10) 




~ .085 






C, = {IS\ Okin 


^^)phys ~ 


0.020(32) 




~ 0.03 






(15 Okin 


\^^)phys 








mQ) 


-E{1S) 


c 


- 




-A 


0.58(4) 


^QCD 


~ 0.36 








1 + 

2 


Okin 


^">. + 

/ phys 


-0.12(8) 




0.1 








Ofiyp 


2 


/ phys 












(15 Okin 


\'^^)phys 










-E{IS) 


Co = 


= A 




-A 


0.40(2) 


^QCD 








Ci = (iPi 


Okin 


1 p \ _ 

^1 phys 


0.36(8) 










(15 Okin 


\'^^)phys 









TABLE V. Predictions for various mass splittings from this work and previous results on 
quenched configurations Q. The experimental values are also included. The results in both cases 
are converted into physical units (MeV) using the scale from rup, given in table 0. The error shown 
is purely statistical, except for Bg—B, where the uncertainty in Kg is shown. The central value for 
this splitting is set using Kg from the K meson. 



/3"/=o = 6.0 /3'^/=2 = 5.6 expt 

B* -B 24(5) 20(12) 45.7(4) 

Bl-Bs 27(3) 26(8) 47(3) 

Bg-B 87(9)(-F20) 98(8)(-12) 90(2) 

5(25) - 5(15) 602(86) 540(120) ~ 580 

Bs{2S) - B,{IS) 559(55) 500(80) 

B** - B 474(32) 520(40) 350-500 

Ab-B 388(68) 560(40) 345(9) 
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TABLE VI. The zeroth order decay amplitude of the PS meson, a^/^f^l^/M^, for all aMo 
and AC;. 



aMo 


0.1385 


0.1393 


0.1401 




0.8 
1.0 
1.2 
1.7 
2.0 
3.0 
3.5 
4.0 
7.0 
10.0 


0.203(3) 
0.214(3) 
0.221(4) 
0.236(4) 
0.243(4) 
0.262(4) 
0.268(4) 
0.274(4) 
0.295(5) 
0.309(6) 


0.190(3) 
0.200(3) 
0.207(4) 
0.220(4) 
0.227(4) 
0.244(5) 
0.251(5) 
0.256(6) 
0.275(6) 
0.288(7) 


0.176(6) 
0.185(6) 
0.191(6) 
0.203(6) 
0.210(7) 
0.227(7) 
0.234(5) 
0.243(8) 
0.270(19) 
278(26) 


0.167(4) 
0.175(5) 
0.181(6) 
0.192(5) 
0.197(6) 
0.212(7) 
0.219(6) 
0.224(8) 
0.239(8) 
249('8) 


TABLE VIL 
all aMo and ki. 


The tree-level correction to the PS current, 




= a^/'fPsVMps, for 


aMo 


0.1385 


0.1393 


0.1401 




0.8 
1.0 

1.7 
2.0 
3.0 
3.5 
4.0 
7.0 
10.0 


-0.0520(10) 
-0.0467(9) 
-0.0422(9) 
-0.0340(7) 
-0.0308(5) 
-0.0234(4) 
-0.0209(4) 
-0.0189(4) 
-0.0120(2) 
-0.0089(2) 


-0.0494(9) 
-0.0443(8) 
-0.0399(8) 
-0.0321(7) 
-0.0291(6) 
-0.0220(5) 
-0.0197(5) 
-0.0178(5) 
-0.0113(3) 
-0.0084(3) 


-0.0462(17) 
-0.0414(15) 
-U.Uo(z(ioj 
-0.0299(10) 
-0.0271(10) 
-0.0207(7) 
-0.0185(4) 
-0.0170(6) 
-0.0112(9) 
-0.0082(8) 


-0.0446(13) 
-0.0398(12) 
-U.Uooo(^iz j 
-0.0286(10) 
-0.0258(10) 
-0.0196(8) 
-0.0175(5) 
-0.0159(6) 
-0.0100(4) 
-0.0074(3) 
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TABLE VIII. The zeroth order decay amplitude of the V meson, a^/2/^°^-/My, for all aMo 
and Ki- 



aMo 


0.1385 


0.1393 


0.1401 




0.8 


0.163(3) 


0.152(3) 


0.140(5) 


0.131(4) 


1.0 


0.176(3) 


0.164(3) 


0.151(5) 


0.142(4) 


1.2 


0.186(3) 


0.174(3) 


0.159(6) 


0.150(4) 


1.7 


0.206(3) 


0.191(4) 


0.177(4) 


0.165(5) 


2.0 


0.215(3) 


0.201(4) 


0.186(4) 


0.173(5) 


3.0 


0.239(4) 


0.223(5) 


0.208(8) 


0.194(6) 


3.5 


0.248(4) 


0.232(5) 


0.216(6) 


0.202(7) 


4.0 


0.255(5) 


0.239(6) 


0.224(8) 


0.209(7) 


7.0 


0.283(6) 


0.264(6) 


0.254(27) 


0.230(10) 


10.0 


0.300(6) 


0.279(7) 


0.272(23) 


0.241(9) 


A "DT TV" 


The tree-level correction to the V current, 


o^^Vy V^y, for all aMo and ki. 


aMo 


0.1385 


0.1393 


0.1401 




0.8 


0.0136(3) 


0.0128(3) 


0.0119(5) 


0.0114(5) 


1.0 


0.0126(3) 


0.0119(3) 


0.0111(5) 


0.0106(4) 




0.0117(3) 


0.0110(3) 


0.0102(4) 


0.0097(3) 


1.7 


0.0098(2) 


0.0092(2) 


0.0086(3) 


0.0082(3) 


2.0 


0.0090(1) 


0.0085(2) 


0.0079(3) 


0.0075(3) 


3.0 


0.0071(1) 


0.0067(2) 


0.0063(3) 


0.0060(2) 


3.5 


0.0064(1) 


0.0061(2) 


0.0057(2) 


0.0054(2) 


4.0 


0.0059(1) 


0.0055(2) 


0.0052(2) 


0.0050(2) 


7.0 


0.0038(1) 


0.0036(1) 


0.0035(4) 


0.0032(1) 


10.0 


0.0029(1) 


0.0027(1) 


0.0027(2) 


0.0024(1) 
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TABLE X. The one-loop correction to the V current, a^/^/^^^/My = a^/^/^^^ y^iWy , for all 
aMo and k;. 



aMo 


0.1385 


0.1393 


0.1401 




0.8 


-0.0164(3) 


-0.0154(3) 


-0.0144(5) 


-0.0136(6) 


1.0 


-0.0144(3) 


-0.0136(3) 


-0.0126(5) 


-0.0119(4) 


1.2 


-0.0129(3) 


-0.0121(3) 


-0.0112(4) 


-0.0107(4) 


1.7 


-0.0103(2) 


-0.0097(2) 


-0.0091(3) 


-0.0086(3) 


2.0 


-0.0094(1) 


-0.0088(2) 


-0.0082(3) 


-0.0078(3) 


3.0 


-0.0072(1) 


-0.0068(2) 


-0.0064(3) 


-0.0061(2) 


3.5 


-0.0065(1) 


-0.0061(2) 


-0.0058(2) 


-0.0055(2) 


4.0 


-0.0059(1) 


-0.0056(2) 


-0.0053(2) 


-0.0050(2) 


7.0 


-0.0038(1) 


-0.0036(1) 


-0.0035(4) 


-0.0032(1) 


10.0 


- 0.0029(1) 


-0.0027(1) 


-0.0027(2) 


-0.0024(1) 


TABLE XI. 


The tree-level and 


renormalised 


PS decay constants 


in lattice units, 


a^^'^ fpsV^PSi for sl\ aMo and = Kc 


aMo has been used for the argument of the logarithms 


appearing in the 


matching coefficients. 








aMo 


tree-level 




q = TT 


q* = 1.0 


0.8 


0.122(3) 




0.121(3) 


0.120(3) 


1.0 


0.135(4) 




0.130(4) 


0.126(3) 


1.2 


0.145(5) 




0.137(4) 


0.132(4) 


1.7 


0.163(5) 




0.150(4) 


0.141(3) 


2.0 


0.172(5) 




0.155(4) 


0.145(4) 


3.0 


0.192(5) 




0.169(5) 


0.154(5) 


3.5 


0.202(5) 




0.175(5) 


0.158(4) 


4.0 


0.208(7) 




0.179(6) 


0.161(5) 


7.0 


0.229(8) 




0.195(6) 


0.174(5) 


10.0 


0.242(9) 




0.206(7) 


0.184(6) 



34 



TABLE XII. The tree-level and renormalised V decay constants in lattice units, a^^"^ fvVMv, 
for all aMf) and ki = Kc- uMq has been used for the argument of the logarithms appearing in the 
matching coefficients. 



aMo 


tree-level 


= TT 


q* = 1.0 


0.8 


0.143(4) 


0.134(4) 


0.128(4) 


1 n 


n 1 A\ 

u. Lo^y^ j 


n 1 A\ 

u. loy j 


u. loi J 


1.2 


0.159(5) 


0.143(4) 


0.134(4) 


1.7 


0.173(5) 


0.152(5) 


0.139(4) 


2.0 


0.181(6) 


0.156(5) 


0.141(4) 


3.0 


0.200(6) 


0.168(5) 


0.148(5) 


3.5 


0.207(7) 


0.173(6) 


0.151(6) 


4.0 


0.214(8) 


0.177(6) 


0.154(6) 


7.0 


0.232(9) 


0.190(8) 


0.164(7) 


10.0 


0.244(9) 


0.199(7) 


0.170(6) 



TABLE XIII. The corrections to the zeroth order pseudoscalar and vector matrix elements ex- 
pressed as a ratio to (^J^^'^'^. The order at which each correction contributes is indicated, aq* = 1.0 
and qMq = 2.0, which is close to the bare b quark mass for = 2.0 GeV. Note that (a and (v 
are set to zero. The statistical errors are less than 1%. 



(J(0)) 

0(1/M) 


apo 
0{a) 


api(j(l)) 
(J(0)) 

0(a/M) 


(J(0)) ' (J{0)^ 

0{a/M) + 0{aa) 


(J(0)) 

0(a/M) 


C4(JW) 
(J(0)) 

0{a/M) 


fpsVMps -13 


-9 


+1 


-7 






fvVMv +4 


-14 


+1 


-6 


-1 


-1 
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TABLE XIV. The pseudoscalar and vector decay constants in MeV, calculated with and with- 
out rcnornialisation and converted to physical units using a^^ = 2.0 and 2.4 GeV. Note that 
aM^ = 2.0 for 0"^ = 2.0 GeV, while aM^ = 1.8 for = 2.4 GeV. Statistical errors, only, are 
shown except for those involving /b^ and fs* , where the error due to the uncertainty in Kg is also 
given. 





a^^ GeV 


tree-level 


q* = IT 


q* = 1.0 


fB 


2.0 


217(6) 


192(5) 


180(5) 




2.4 


269(8) 


244(8) 


231(8) 




2.0 


247(5)(-5) 


222(2) (-5) 


208(2) (-5) 




2.4 


310(6) (-6) 


284(4) (-6) 


266(4) (-6) 


fs* 


2.0 


223(6) 


189(6) 


173(5) 




2.4 


280(8) 


248(8) 


224(6) 


fB* 


2.0 


256(4) (-3) 


224(4) (-4) 


202(4) (-4) 




2.4 


323(6) (-6) 


285(5) (-7) 


261 (5) (-5) 


fBs/fB 




1.14(2)(-2) 


1.14(2)(-2) 


1.14(2)(-2) 


/b*/ /b* 




1.14(2)(-2) 


1.14(2)(-2) 


1.14(2)(-2) 



TABLE XV. The dependence of the 0(1/M) slope of a^^^f^/M^ on the light quark mass in 
lattice units. 



Hi 


aCi/Co 


0.1385 


1.7(5) 


0.1393 


1.0(4) 


0.1401 


1.1(3) 




1.1(3) 



TABLE XVL The slope of a^^"^ fy/Mps in lattice units at tree-level and when renormalisation 
is included. 





a3/2Co 


a5/2Ci 


aCi/Co 


tree-level 


0.27(1) 


-0.3(1) 


-1.1(3) 


q = TT 


0.192(8) 


-0.09(2) 


-0.5(1) 


q* = 1.0 


0.150(6) 


-0.00(2) 


-0.0(1) 
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TABLE XVII. The 0(1/M) slope found for various combinations of the tree-level PS and V 
decay constants in lattice units, ki = Kc- 



Quantity 




aCi/Co 


Results 


JPSVMps 
fVM 

fpsVMps 




Gkin 


-1.1(3) 
-0.9(3) 

0.36(3) 

-0.14(4) 


Quantity 




aCo 


Results 


ffl^Mps 




Gcorr 


-0.62(1) 


TABLE XVIII. A comparison of the decay constants of the B and Bs mesons from this work 
and those on quenched configurations at j3'^f^^ = 6.0 The results have been converted into 
physical units (MeV) using from nip. The errors are purely statistical, with the exception of 
fSs/fB where the second error comes from the uncertainty in Kg. Note that the 0(1/M^) current 
corrections included in reference Q have been omitted from the quenched results in order to better 
compare with the nj = 2 results. 






/3"/=o = 6.0 


/3"/=2 = 5.6 


fB 

fSs 1 fs 


q* = 1.0 
q* =7: 
tree-level 
q* = 1.0 
q* =7r 
tree-level 


145(10) 
154(10) 
176(11) 
175(6) (+7) 
186(6)(+7) 
205(6) (+7) 
1.20(4)(+4) 


180(5) 
192(5) 
217(6) 
208(2) (-4) 
222(2) (-4) 
247(2) (-4) 
1.14(2)(-2) 
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TABLE XIX. The operators and corresponding quantum numbers used for the computation of 
the spectrum of the heavy-light mesons. 



2S+1 



Lj (J 



PC\ 



'So (0-+) 

(1--) 

'Pi (1+-) 
'Po (0++) 
'Pi (1++) 

(2++) 



q^iAiaj - Ajai)Q 
q^{Aiai - Ajaj)Q 
q^AiUj + Ajai)Q {i^j) 



TABLE XX. The ground state pseudoscalar meson energies, aEQ-So), extracted from Uexp = 2 
fits to C{1, 1) and C{1, 2) correlators. The mass shift needed to convert the simulation energies to 
the meson mass is also given. 



aMo 


A 


0.1385 


0.1393 


0.1401 




0.8 


0.89(2) 


0.420(4) 


0.395(4) 


0.367(7) 


0.347(7) 


LO 


1.09(3) 


0.465(4) 


0.440(5) 


0.412(8) 


0.394(8) 


1.2 


1.27(2) 


0.492(4) 


0.467(5) 


0.439(8) 


0.421(9) 


1.7 


1.76(2) 


0.524(5) 


0.500(5) 


0.473(8) 


0.453(9) 


2.0 


2.07(2) 


0.534(5) 


0.509(5) 


0.483(8) 


0.463(10) 


3.0 


2.58(2) 


0.548(5) 


0.524(6) 


0.499(9) 


0.478(7) 


3.5 


3.10(3) 


0.551(5) 


0.528(6) 


0.503(9) 


0.483(9) 


4.0 


3.64(3) 


0.553(5) 


0.530(7) 


0.507(9) 


0.486(9) 


7.0 


6.61(12) 


0.558(7) 


0.536(7) 


0.515(10) 


0.492(10) 


10.0 


9.43(22) 


0.560(8) 


0.536(8) 


0.515(11) 


0.491(12) 
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TABLE XXI. The first excited state pseudoscalar meson energies, aE{2^So), extracted from 
Uexp = 2 fits to C{1, 1) and C{1,2) correlators. 



aMo 


0.1385 


0.1393 


0.1401 




0.8 


0.70(5) 


0.68(5) 


0.68(6) 


0.68(8) 


1.0 


0.73(4) 


0.72(5) 


0.71(5) 


0.70(7) 


1.2 


U.75(oj 


U.74(4j 


0.7z(5j 


0.71(6j 


1.7 


0.77(3) 


0.76(4) 


0.74(5) 


0.73(6) 


2.0 


0.78(3) 


0.77(4) 


0.74(5) 


0.73(5) 


3.0 


0.78(3) 


0.77(3) 


0.74(4) 


0.73(5) 


3.5 


0.78(3) 


0.76(3) 


0.74(5) 


0.73(5) 


4.0 


0.77(3) 


0.76(3) 


0.74(4) 


0.73(5) 


7.0 


0.76(3) 


0.75(3) 


0.75(4) 


0.73(5) 


10.0 


0.74(3) 


0.73(3) 


0.73(5) 


0.72(5) 


TABLE XXII. The ground state vector meson energies, 


aE[ Si), extracted from 




to C{1, 1) 


and C(Z,2) correlators. 








aMo 


0.1385 


0.1393 


0.1401 




0.8 


0.448(4) 


0.422(5) 


0.391(8) 


0.372(9) 


1.0 


0.490(4) 


0.463(5) 


0.432(9) 


0.415(9) 


1.2 


0.513(4) 


0.487(6) 






1.7 


0.541(4) 


0.515(6) 


0.485(9) 


0.470(9) 


2.0 


0.548(4) 


0.523(5) 


0.493(9) 


0.477(8) 


3.0 


0.558(4) 


0.534(5) 


0.507(9) 


0.488(8) 


3.5 


0.560(5) 


0.536(5) 


0.510(8) 


0.490(7) 


4.0 


0.561(5) 


0.537(5) 


0.512(8) 


0.492(8) 


7.0 


0.564(6) 


0.540(7) 


0.518(11) 


0.494(9) 


10.0 


0.565(10) 


0.541(9) 


0.519(12) 


0.495(13) 
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TABLE XXIII. The first excited state vector meson energies, aE{2^Si), extracted from 
^exp = 2 fits to C{1, 1) and C{1,2) correlators. 



aMo 


0.1385 


0.1393 


0.1401 




0.8 


0.75(5) 


0.74(6) 


0.71(6) 


0.69(7) 


1.0 


0.78(5) 


0.77(6) 


0.74(6) 


0.72(7) 


1.2 


n on / /I ^ 

U.oU(4j 


0.79(5; 


0.75(5j 


0.7z(6j 


1.7 


0.81(4) 


0.81(4) 


0.76(5) 


0.75(6) 


2.0 


0.82(3) 


0.81(4) 


0.76(5) 


0.73(6) 


3.0 


0.81(3) 


0.80(4) 


0.76(5) 


0.74(6) 


3.5 


0.80(3) 


0.79(4) 


0.76(4) 


0.74(5) 


4.0 


0.80(3) 


0.79(4) 


0.76(4) 


0.75(5) 


7.0 


0.77(3) 


0.76(3) 


0.75(4) 


0.75(5) 


10.0 


0.74(3) 


0.74(4) 


0.74(4) 


0.73(5) 


TABLE XXIV. The ground state enerj 


of the ^Pi state, aE(^Pi), for various aMo and all ki. 


aMo 


0.1385 


0.1393 


0.1401 


He 


1.0 


0.79(1) 


0.77(1) 


0.75(1) 


0.74(1) 


2.0 


0.78(1) 


0.76(1) 


0.75(1) 


0.74(1) 


4.0 


0.76(1) 


0.75(1) 


0.74(1) 


0.71(1) 


TABLE XXV. The ground state energy of the Aq 


, aE^Ag), extracted from Uexp 


= 1 fits to 


C{1,1) 


for ah aMo and k;. 








aMo 


0.1385 


0.1393 


0.1401 




0.8 


0.783(9) 


0.727(13) 


0.669(17) 


0.621(20) 


1.0 


0.826(10) 


0.773(13) 


0.716(17) 


0.671(18) 


1.2 


0.851(10) 


0.799(12) 


0.744(17) 


0.701(17) 


1.7 


0.880(9) 


0.832(11) 


0.779(17) 


0.739(18) 


2.0 


0.889(9) 


0.841(10) 


0.789(19) 


0.751(18) 


3.0 


0.902(9) 


0.858(10) 


0.804(17) 


0.772(18) 


3.5 


0.905(9) 


0.862(10) 


0.808(18) 


0.778(18) 


4.0 


0.907(9) 


0.865(12) 


0.811(19) 


0.782(17) 


7.0 


0.893(17) 


0.873(13) 


0.821(18) 


0.801(28) 


10.0 


0.886(23) 


0.878(12) 


0.828(24) 


0.819(26) 
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FIG. 1. The lower lying B spectrum and A;,. The open circles denote our results. The dashed 
lines denote the upper and lower bounds on the experimental results. The dotted lines indicate 
preliminary experimental signals. For the P states the dashed and dotted lines represent the B^*^Tr 
and the narrow Btttt resonance respectively. The errors shown are purely statistical. The filled 
circles represent preliminary results from a quenched simulation at (3''^^^^ = 6.0 Q. All results are 
converted to physical units using the scale from irip, given in table ||. 
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FIG. 2. The theoretical prediction for the structure of the B spectrum [27|. ji and J denote 
the total spin of the light quark and the heavy-light meson respectively. 
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0.0 I I I I I ' ' ' I I I I I ' I ' ' ' I I I ' I ' 
0.0 0.2 0.4 0.6 0.8 1.0 
l/(aMps) 

FIG. 3. The dependence of the mass splittings of the ground and first excited state PS and V 
mesons with the inverse pseudoscalar meson mass in lattice units, m = Kc for the E{2S)—E{1S) 
splitting, while ki = 0.1385 for the IS and 2S hyperfine splittings. 
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FIG. 4. The dependence of the E{^So)^^-E{^So)^^, E{Aq)-E{1S) and E{1P) - E{1S) spht- 
tings on l/{aMps) in lattice units, where m = Kc- 
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0.35 



1.20 



S: 



"a 0.20 — 



0.10 

0.0 0.2 0.4 0.6 O.i 

l/(aMps) 



(b) 



0.0 0.2 0.4 0.6 O.f 

l/(aMps) 



FIG. 5. (a) shows the tree-level a^^'^f^/Mps for the three values of ki = 0.1385 (crosses), 
0.1393 (octagons) and 0.1401 (diamonds), and extrapolated to the chiral limit (squares) as a 
function of l/{aMps). (b) shows the heavy quark mass dependence of the ratio of {f\fMps)Ks to 
{f^/Mps)Kci where Kg is fixed using the K meson. 
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FIG. 6. The pseudoscalar decay constant in lattice units fully consistent to 0{a/M) as a func- 
tion of l/{aMps). The results shown as circles are obtained using aq* = 1.0 for the renormalisation 
factors, while the squares use aq* = vr; the trcc-lcvcl results are also shown as diamonds. aMg = 2.1, 
for = 2.0 GeV, has been used for the argument of the logarithms appearing in the matching 
coefficients, for all gMq. The results are for ki = Kc- 
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FIG. 7. The ratio of individual tree-level and 1-loop current corrections to the zeroth order 
PS decay constant, /("^v/M^, as a function of l/(aMo). = (Jl^) / (Jl'')- i^l = i^c and 

aq* = 1.0. aMo = qMq has been used for the argument of the logarithms appearing in the 
matching coefficients. 
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FIG. 8. The individual tree-level and 1-loop current corrections to the zeroth order V decay 
constant, fi^^VM^ as a function of l/(aMo). R^^^ = (4*'*) / (4°^)- = '^c and aq* = 1.0. 
aMo = aM^ has been used for the argument of the logarithms appearing in the matching coeffi- 
cients. 
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FIG. 9. Heavy quark dependence of various combinations of the tree-level PS and V decay 
constants in lattice units for ki = k^. (a) presents fyJMps (diamonds), f^^^ yJMps (octagons) and 
the spin-average of the PS and V decay constants (squares), (b) shows the ratio of the PS to the V 
tree-level decay constant with (diamonds) and without (squares) the 0(1/M) current corrections, 
(c) gives the contribution to the slope of the decay constants from the 0(1/M) current correction. 
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FIG. 10. The renormalised pseudoscalar decay constant as a function of l/Mpg in physical 
units for /J"/^*^ = 6.0 (diamonds 1/M^) and = 5.6 (squares) for tree-level and 1-loop for 

q* = TT and q* = 1.0. aMg has been used for the argument of the logarithms appearing in the 
matching coefficients. 
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FIG. 11. The ground and first excited state energies extracted as a function of tmin from a 
vector fit to tire C{1, 1) and C{1,2) pseudoscalar correlators for uMq = 1.0 and ki = 0.1385; tmax 
is fixed to 20. E{^So) is shown as squares (crosses) for a one (two) exponential fit and E{2^So) as 
circles for the two exponential fit. The quality of fit parameter, Q, is also shown as a dotted and 
solid line for a one and two exponential fit respectively, where Q > 0.1 defines a 'good' fit. 
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FIG. 12. The chiral extrapolation of aE{^So) and a^/"^ f^^^ y/M^ in 1/k for aMo = 1.0. 
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FIG. 13. The ground state energy of the ^Pi state extracted from a n^xp = 1 fit to C{1,2) as 
a function of tfnin for clMq = 1.0 and ki = 0.1385; tmaa; ~ 

10. The corresponding values of Q are 

also shown. 
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FIG. 14. The ground state energy of the Aq extracted from a Ugxp = 1 (crosses) and 
Uexp = 2 (circles) fit to the C{1, 1) correlator as a function of tmin for aMo = 1.0 and ki = 0.1385; 
tmax = 20. The corresponding values of Q are also shown as a solid (dotted) line for a Uexp = 1 (2) 
fit. 
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FIG. 15. The ground state energy extracted from a n^xp = 1 fit to the pseudoscalar C{1, 1) and 
C(l, 1) meson correlators as a function of tmin for oMq = 1-0 and = 0.1385; tmax = 20. The 
corresponding values of Q are also shown. 
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FIG. 16. The ratio C{lj(i),l)/C{l,l) for the PS meson as a function of tmin for clMq = 1.0 
and Ki = 0.1385; tmax = 20. 
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